Mon. Not. R. Astron. Soc. 000 , 1-14 (2008) 


Printed 29 April 2008 (MN lATftjX style file v2.2) 


Is there evidence for X-ray emitting plasma very close to 
the photospheres of O stars? 

Maurice A. Leutenegger 1,2,3,4 

1 Laboratory for High Energy Astrophysics, Code 662, NASA/Goddard Space Flight Center, Greenbelt, MD 20771 

2 Department of Physics and Astronomy, Swarthmore College, Swarthmore, PA 19081 

3 Department of Physics and Columbia Astrophysics Laboratory, Columbia University, 550 West 120th Street, New York, NY 10027 
4 NASA Postdoctoral Fellow 

Submitted 2008 Xxxxx XX 


ABSTRACT 

I describe a procedure for accurately assessing parameter confidence limits for 
low signal-to-noise X-ray emission line complexes in which the blending of Doppler- 
broadened lines is significant. I apply this methodology to measure the f/i ratios 
of high-Z He-like ions in the X-ray spectra of several O stars. I find that the data in 
question actually provide no significant constraint on the location of the X-ray emitting 
plasma. This is not surprising, as the data are of poor statistical quality. The data thus 
do not require a two-component model for the origin of O-star X-ray emission, and do 
not present a near-star high-ion problem, as suggested in Waldron & Cassinelli (2007). 
This is consistent with earlier results, using lower-Z He-like ions, constraining the bulk 
of the X-ray emitting plasma to be distributed throughout the wind and not near the 
photosphere. It is also consistent with the consensus finding that emission lines from 
high-Z ions in the spectra of many O stars have significant Doppler broadening, which 
is not expected for lines forming very close to the photosphere. 
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1 INTRODUCTION 

The first high resolution X-ray spectra of O stars obtained 
by the diffraction grating spectrometers onboard XMM- 
Newton and Chandra showed He-like triplet emission with 
strong intercombination lines ( i ) and weak or absent forbid- 
den lines (/) . This was initially surprising; many researchers 
were familiar with the idea that this was a signature of a 
high-density plasma, and the densities implied would have 
been difficult to understand if they were present in the winds 
of O stars, where the X-rays are thought to originate. In their 
investigation of the RGS spectrum of £ Puppis, Kahn et al. 
(2001) showed that the observed line ratios were due to the 
strong ultraviolet radiation from the photosphere, which ex- 
cites electrons from the metastable ls2s 3 Si state (the upper 
level of the forbidden line) to the ls2p 3 Pj states (the up- 
per levels of the intercombination lines). The dependence of 
the f/i ratio on the UV field was in fact derived in the first 
works on this subject (Gabriel & Jordan 1969; Blumenthal 
et al. 1972), but, in the absence of an obvious context for its 
importance, it was less widely appreciated than the density 
dependence. The dependence of the observed ratio on the 
strength of the UV field allowed Kahn et al. (2001) to use 


the ratios to infer rough plasma locations. They found that 
X-rays were indeed originating from the wind of C, Pup. 

Waldron & Cassinelli (2001) subsequently used the for- 
malism of Gabriel & Jordan (1969) and Blumenthal et al. 
(1972) to derive relatively precise radii of formation from the 
observed f/i ratios for Orionis, under the assumption of a 
single formation radius for a given ion. They found that most 
of the X-ray emitting plasma was located at radii consistent 
with a wind-shock origin. However, one of the more perplex- 
ing results of this investigation was the inference from the 
f/i ratio of Si XIII of the presence of hot X-ray emitting 
plasma very close to the photosphere of (j Ori (at a height 
of less than 0.1 stellar radii). Cassinelli et al. (2001) made 
a similar inference based on the f/i ratio of S XV observed 
in the X-ray spectrum of () Pup, and Waldron et al. (2004) 
again found this behavior for S XV and Ar XVII in the X-ray 
spectrum of Cyg OB2 8A. Again, in each of these cases, the 
bulk of the X-ray emission was found to originate from loca- 
tions consistent with shocks distributed in the wind, while 
only the highest-Z He-like ions were inferred to be forming 
near the photosphere. 

Although the initial prediction of X-ray emission from 
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O stars based on the observed super-ion states was couched 
in the context of a base coronal model (Cassinclli & Olson 
1979), both the UV observations of super-ions and the subse- 
quent detection of X-rays from O stars (Harnden et al. 1979; 
Seward et al. 1979) indicate nothing about the mechanism 
of X-ray production. Moderate resolution X-ray spectra ob- 
tained with the Einstein Solid State Spectrometer showed 
that at least the soft X-ray emission must be distributed 
in the winds of O stars and cannot be located in a corona 
or other confined region near the photosphere (Cassinclli & 
Swank 1983). Around the same time, it was proposed by 
Lucy & White (1980) that the instability of the radiative 
driving mechanism of O star winds leads to the formation 
of shocks which are capable of producing X-ray emission 
throughout the winds of O stars. A large body of theoret- 
ical work has further developed this idea (e.g., Lucy 1982; 
Macfarlane et al. 1993; Cooper 1994; Cooper & Owocki 1994; 
Feldmeier et al. 1997,?; Owocki & Cohen 1999). 

The spectra of most O stars observed with the diffrac- 
tion grating spectrometers onboard XMM - Newton and 
Chandra are consistent with the wind-shock paradigm for X- 
ray plasma formation. The strongest single piece of evidence 
for this is the dramatic velocity broadening of the emission 
lines (including the lines of the He- like triplets); they typ- 
ically have HWHM of order half the wind terminal veloc- 
ity, or about lOOOkms -1 (e.g., Kahn et al. 2001; Cassinclli 
et al. 2001; Kramer et al. 2003; Cohen et al. 2006). This is 
in contrast with coronal sources, which never have velocity 
broadening or shifts that are detectable at the resolution of 
the current generation of instruments (<il00 kms -1 ), except 
for orbital motion in a few close binaries, which is barely de- 
tectable. Although there are still unsettled questions in the 
interpretation of O star line profiles, especially regarding 
the relative importance of mass-loss rate reduction, poros- 
ity, and resonance scattering (e.g., Owocki & Cohen 2001; 
Ignace & Cayley 2002; Kramer et al. 2003; Oskinova et al. 
2004; Cohen et al. 2006; Owocki & Cohen 2006; Oskinova 
et al. 2006; Leutenegger et al. 2007), there is no doubt that 
the non-negligible line widths should be interpreted in the 
context of some kind of wind-shock model. 

A few anomalous O and early B stars are clearly not 
consistent with the instability wind-shock paradigm for X- 
ray emission; in particular, 8 1 Orionis C (Gagne et al. 2005) 
and r Scorpii (Cohen et al. 2003). They are overluminous 
and have higher temperature plasma compared to other O- 
star X-ray sources, and their emission lines are relatively 
narrow, although they are resolved. In addition, each of these 
two stars has a relatively strong magnetic field (Donati et al. 
2002, 2006). ud-Doula & Owocki (2002) have performed nu- 
merical magnetohydrodynamic simulations of the wind of 
8 1 Ori C, and they find that the magnetic field channels 
part of the wind and causes it to collide at the magnetic 
equator with high relative velocity, thus explaining the ob- 
served high X-ray temperature and luminosity. It is likely 
that a mechanism involving the interaction of the wind and 
magnetic field of r Sco is also responsible for its high X-ray 
luminosity and high temperature. 

The fact that these are relatively unevolved stars with 
non-negligible large-scale magnetic fields has led to a pic- 
ture in which young OB stars may have fossil magnetic 
fields which channel their winds and create relatively high- 
luminosity, hard X-ray emission, while evolved OB stars no 


longer have strong magnetic fields and generate soft X-rays 
through the line-driven-instability wind-shock mechanism. 
The observational evidence and theoretical work supporting 
this picture is reviewed in Cohen (2008). 

The presence of any X-ray emitting material close to the 
stellar photosphere of a “normal” O star would be difficult 
to explain in the context of a wind-shock model. Numeri- 
cal hydrodynamic simulations of instabilities in radiatively 
driven stellar winds indicate that several tenths of a stel- 
lar radius are required to allow the instabilities to develop 
(e.g., Owocki et al. 1988) if the line-driving instability is 
self-excited; the presence of perturbations at the base of the 
wind might alleviate the problem to some degree, but gener- 
ating X-rays very close to the photosphere (less than a tenth 
of a stellar radius) would likely still be a problem (Feldmeier 
1995). However, note that numerical simulations give no in- 
dication that the characteristic shock velocity scales with 
the local flow velocity, as suggested by Waldron & Cassinelli 
(2001), Cassinelli et al. (2001), and Waldron & Cassinelli 
(2007). On the contrary, these simulations typically show 
reverse shocks caused by rapidly accelerating flows colliding 
with clumps that follow the trajectory of the mean wind 
velocity law, with large shock velocities even for the shocks 
occuring closest to the star (Owocki et al. 1988; Feldmeier 
1995; Runacres & Owocki 2002). 

The inference of X-ray emitting plasma forming close to 
the photospheres of the O stars ( Ori (Waldron & Cassinelli 
2001), £ Pup (Cassinelli et al. 2001), and Cyg OB2 8A (Wal- 
dron et al. 2004) led Mullan & Waldron (2006) to suggest 
the possibility of a two-component origin for their X-ray 
emission, with the bulk of the X-ray emission arising from 
shocks distributed throughout a wind flowing from the equa- 
tor, but with the highest temperature plasma originating in 
polar coronal loops very close to the photosphere. The high- 
ion f /i ratios which led to the inference of X-ray emitting 
plasma close to the stellar photospheres of those stars are 
the strongest piece of observational evidence presented in 
support of this scenario. 

In Leutenegger et al. (2006), my collaborators and I at- 
tempted to reproduce some of these early results and found 
that the f /i ratios of Si XIII in the X-ray spectrum of £ Ori 
and S XV in the spectrum of ( Pup did not require the pres- 
ence of hot plasma very close to the photosphere, but were 
consistent with formation radii several tenths of a stellar ra- 
dius above the photosphere at the 68% confidence level. We 
identified two main sources of error in those measurements 
of characteristic radii of formation by Waldron & Cassinelli 
(2001) and Cassinelli et al. (2001) that indicated the exis- 
tence of X-ray emitting plasma close to the photosphere: the 
measurements of the line ratios themselves, and the model 
photospheric UV flux shortward of the Lyman break, which 
is relevant for Si XIII and higher-Z He-like ions. 

Recently, Waldron & Cassinclli (2007) (hereafter 
WC07) have published measurements of line ratios as well 
as velocity shifts and widths for a large sample of O stars 
observed with the Chandra High Energy Transmission Grat- 
ing Spectrometer (HETGS). Using f/i line ratios, they have 
inferred radii of formation near the photosphere for high- 
Z lielium-like ions observed in the spectra of a number of 
stars (typically S XV or Si XIII), and they have character- 
ized their observations as constituting a near-star high-ion 
problem. The fact that this behavior is observed for a large 


© 2008 RAS, MNRAS 000 , 1-14 



High-Z-ion f/i ratios in O stars 3 


sample of stars lends greater weight to the claim of a uni- 
versal problem than observations in only three stars. If their 
observations were confirmed, they would present a signifi- 
cant challenge to the wind-shock mechanism for producing 
X-rays from O stars. 

WC07 have published an erratum in which they have 
made significant corrections to the line ratio measurements 
and their confidence intervals. In particular, due to a bug in 
their line-fitting code, the confidence intervals published in 
their original work were significantly underestimated. Fur- 
thermore, they revised their fit parameter confidence inter- 
vals from 90% to 68%. However, they report in the erratum 
that this does not affect their general conclusion that there 
is a near-star high-ion problem. I discuss their measurements 
in more detail in § 3. 

In this paper I attempt to reproduce the measurements 
which lead WC07 to infer radii of formation near the photo- 
sphere for high-Z ions. In all cases that I consider, I find that 
the data provide no significant constraint on the radii of for- 
mation. However, this does not invalidate measurements for 
lower-Z ions, including those of WC07, which have generally 
found that those ions are not located near the photosphere. 

This paper is organized as follows. In § 2 I give a brief 
review of the atomic physics of He-like triplets and the IZ = 
f /i ratio. In § 3 I summarize the f /i ratio measurements of 
WC07. In § 4 I discuss the reduction of the data, present my 
measurements of f/i, and discuss possible explanations for 
the discrepancies between my results and those of WC07. In 
§51 discuss my results, and in § 6 I give a brief summary. 


2 RELATION BETWEEN PLASMA 

FORMATION RADIUS AND F/I RATIO 

In this section I give a brief review of the behavior of the 
7 Z = f/i ratio of a He-like ion in the presence of a strong 
external UV held, as in the wind of an O-star. For a more 
complete discussion, see e.g. Leutenegger et al. (2006). I use 
the notation of that paper throughout. In particular, I adopt 
the use of script 7 Z to denote the f/i ratio, and italic R to 
denote quantities comparable to a stellar radius. 

The diagnostic utility of this ratio originates in the fact 
that the upper level of the forbidden line (/), 1s2s 3 Si, is 
metastable. If a process excites electrons in this state to the 
upper levels of the intercombination lines (*), ls2p 3 Pi, 2 , the 
f/i ratio is altered. The criterion for significant alteration 
of the line ratio is that the excitation rate from 1s2s 3 S'i to 
ls2p 3 Pi , 2 should be comparable to the decay rate of ls2s 3 Si 
to ground. 

There are two processes that can excite electrons in this 
way: electron impact excitation, and photoexcitation. The 
former becomes important at high density, the latter when 
the UV field intensity is strong. 

The 1Z ratio is an especially powerful diagnostic because 
the decay rate to ground of the upper level of the forbidden 
line is a strong function of Z , so that even though the range 
of Z typically observed only spans factor of two, the range 
of densities or UV held strengths probed by the observed 1Z 
ratios spans a few orders of magnitude. 

Gabriel & Jordan (1969) and Blumenthal et al. (1972) 
solved the rate equations to show that 


7 Z{<j>, n e ,T) = IZq(T) 


1 

1 + </>/</>c. + n e /n c (T) ’ 


(1) 


where 1Z = f/i; IZo is the ratio in the limit of low UV 
held strength and low density; <j> is the photoexcitation rate 
from 1 s2s 3 Si to ls2p 3 Pj; n e is the electron density; and (j> c 
and n c are the “critical” photoexcitation rate and density, 
respectively, where 7 Z drops to half of IZo- ij> c , n c , and IZo are 
functions of atomic parameters. n c and IZo depend weakly 
on temperature. 

In the case of an ion in a O-star wind, the value of 
4> is proportional to the geometrical dilution W(r) = |{1 — 
[1 — (Tth/r) 2 ] 1 / 2 }, and the density is negligible in comparison 
with the critical density, so we may write 


U n °l+2 PW{r) 


(2) 


where P is a constant that depends only on atomic physics 
and the photospheric UV hux. By assuming that emission 
from a given ion comes from a single radius, one may invert 
this expression for an observed value of TZ to find a char- 
acteristic radius of formation; this is the approach taken by 
WC07. Other approaches based on more realistic assump- 
tions may be taken to extract information regarding the 
radial location of the plasma (e.g., assuming a radially dis- 
tributed plasma, as in Leutenegger et al. 2006). 

The theory developed in Gabriel & Jordan (1969) and 
Blumenthal et al. (1972) is adequate for most practical ap- 
plications. However, it does not include the effect of blended 
dielectronic recombination (DR) satellite lines from Li-like 
species, which can be important for high-Z He-like ions. 
Given that the temperature distribution of X-ray emitting 
plasmas observed on O stars is skewed to lower tempera- 
tures, this effect is likely significant for Si XIII and higher-Z 
ions. 

As the f/i ratios measured in this paper are poorly 
constrained, I do not use the the formalism described in 
this section to calculate allowed ranges of plasma forma- 
tion radii, and it is sufficient to note that the satellite lines 
will lower the effective observed value of IZo for high-Z ions. 
This effect has been included in the calculations of Porquet 
et al. (2001), but unfortunately the highest.-Z ion considered 
in that work is Si XIII. Thus, for Si XIII, I adopt 77-0 (Tm) 
from Porquet et al. (2001); for S XV, I adopt the value of 
Blumenthal et al. (1972); and for Ar XVII, I interpolate the 
Blumenthal et al. (1972) values for S XV and Ca XIX. Be- 
cause they do not include DR satellites, the IZo (T m ) values 
of Blumenthal et al. (1972) for high-Z ions are probably 
overestimated by at least 10%, but given the low precision 
of the measurements considered in this paper, that is an ad- 
equate level of accuracy. The values of IZo {Tm) adopted are 
summarized in Table 1. 


3 SUMMARY OF F/I RATIO 

MEASUREMENTS IN IN WC07 

In § 3.1 I present a summary of the measurements of WC07 
implying the formation of X-ray emitting plasma close to 
the photospheres of some O stars. In § 3.2 I discuss their 
fitting method and error analysis procedures. 
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Table 1. Adopted values of TZq (T m) 


Ion 

77o (7m) 

Si XIII 

2.3 “ 

S XV 

•O 

O 

CM 

Ar XVII 

1.7 c 


“Porquet et al. (2001) 

^Blumenthal et al. (1972) 

c Extrapolated from Blumenthal et al. (1972) 

3.1 Measurements 

In Table 2 I have collected all f /i ratios reported in WC07 
which led them to infer plasma formation radii close to the 
stellar surface. For purposes of this paper, I have chosen 
the criterion that the reported upper limit to the plasma 
location is < 1.7 R* at 68% confidence. The results given 
are from the erratum of WC07 and not from the original 
publication. 

This does not imply that X-ray emitting plasma as dis- 
tant as 1.7 R* would necessarily be problematic if it were de- 
tected. As discussed in the introduction, numerical modeling 
of the line-driving instability indicates that a wind seeded 
with perturbations could develop strong shocks significantly 
lower than 1.5 R*. The choice to include data constraining 
X-ray emitting plasma to be closer than 1.7 R* is made to 
ensure that all data indicating a possible near-star high-ion 
problem are examined. 

I only show results for Si XIII and S XV, and not for 
lower-Z ions. I also give the results of Waldron et al. (2004) 
for Ar XVII in the spectrum of Cyg OB2 8A, which indicate 
that this ion is also formed close to the photosphere for this 
star. 

WC07 report separate measurements using the Medium 
Energy Grating (MEG) and High Energy Grating (HEG), 
rather than fitting MEG and HEG data simultaneously. The 
MEG and HEG are similar transmission gratings with a fac- 
tor of two difference in dispersion; all HETGS observations 
give simultaneous MEG and HEG spectra. The HEG and 
MEG results are both given for every ion considered, when 
available. In cases where measurements from one grating im- 
ply formation close to the photosphere while those from the 
other do not, I give both for comparison, although there is 
no reason to suppose that MEG and HEG data should give 
results whose differences are statistically significant. 

Of the measurements listed in Table 2, three show con- 
sistent results for both the MEG and HEG measurements, 
two show marginally consistent results, two were only made 
with the MEG and not with the HEG, and one shows in- 
consistent results between the MEG and HEG, given the 
reported 68% confidence intervals. 

In Table 2 I also report the total number of counts in the 
fit range for each complex for both the HEG and the MEG. 
Although some of the complexes have a significant number 
of counts, some of them have so few counts that it is hard to 
see how one could extract any meaningful information from 
them other than the overall normalization. For example, the 
S XV triplet of <( Ori has a total of only 24 counts (MEG 
and HEG combined), of which about half are expected to 
be distributed between the forbidden and intercombination 
lines, with the remainder being in the resonance line. 


Of the 17 stars considered in WC07, ten show no evi- 
dence requiring X-ray emitting plasma close to the photo- 
sphere (Cyg OB2 9, 5 Ori, e Ori, HD 150136, i Ori, (5 Cru, 
9 Sgr, 15 Mon, £ Oph, and r Sco), while seven stars have 
some evidence requiring this. In each of these seven cases, 
only the highest-Z He-like species observed shows any evi- 
dence requiring X-ray emitting plasma close to the photo- 
sphere, except in the case of Cyg OB2 8A, for which both 
S XV and Ar XVII require this, and a Ori, for which both 
Mg XI and Si XIII require this. 

Three He-like complexes that were claimed to require X- 
ray emitting plasma close to the photosphere in the original 
version of WC07 (S XV in the spectrum of ( Oph, and Si XIII 
in the spectra of t Ori and ( Ori) are no longer claimed to 
require this in the erratum. This is particularly significant 
in the case of ( Ori. The formation radius of Si XIII inferred 
by Waldron & Cassinclli (2001) for this star was the ori- 
gin for the first claim of a near-star high-ion problem, and 
led Waldron & Cassinclli (2001) to infer that some of the 
X-ray emitting plasma was located close enough to the pho- 
tosphere that high-density effects might become important. 
This claim is implicitly rescinded by the new measurements 
of WC07 published in their erratum: the HEG data now 
give no constraint, and the MEG data give a lower limit of 
R > 1.8R * . 


3.2 Fitting method and error analysis 


In this section I give a brief summary of the methodology of 
WC07 based on their original publication, the erratum, and 
private communications from W. Waldron. 

The data in the spectral region of each He-like complex 
are modeled by three Gaussian emission lines plus a contin- 
uum component. Each Gaussian has its own normalization, 
but they all share the same centroid Doppler shift, and the 
same Doppler width. 

The data are fit by finding the point in model parameter 
space that minimizes the statistic 


2 \ ' (Ni — Mi ) 2 

X a 2 


( 3 ) 


where i are the indices of spectral bins, Ni are the observed 
counts in each bin, M, are the predicted model counts, and 
<t; are the variances associated with the observed counts. 
The variances are computed using 

at = max[1.87, VNi] , (4) 


This function has similar behavior to the Gehrels variance 
function (Gehrels 1986), 

at = 1 + V0. 75 + Ni , (5) 

which can be used in XSPEC or other similar spectral fitting 
programs. 

It is not clear how parameter confidence intervals were 
computed in WC07. The standard way to estimate model 
parameter errors is to vary the parameters and compare 
the change in fit statistic due to the variation (e.g. Lamp- 
ton et al. 1976; Press et al. 2007, “Confidence Limits on 
Estimated Model Parameters”, Chapter 15.6 in Numerical 
Recipes). 

For example, to find the confidence intervals for 77, one 
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Table 2. Measured 7 Z ratios reported in WC07 that imply formation radii close to 
the photosphere 


MEG HEG 


Star 

Ion 

77° 


R b 

N c 

n a 


R b 

N c 

Cyg OB2 8A 

Ar: 

xvn d 

< 

0. 

16 


1.0 

64 

0.43 

± 

0.37 


1.0 

53 

Cyg OB2 8A 

s 

XV 

0.74 

± 

0.34 

< 

1.14 

140 

0.87 

± 

0.37 

< 

1.26 

90 

C Pup 

s 

XV 

0.60 

± 

0.62 

< 

1.48 

85 

0.79 

± 

0.65 

< 

1.84 

46 

e 1 Ori C 

s 

XV 

1.32 

± 

0.66 

< 

4.88 

495 

0.70 

± 

0.25 

< 

1.04 

365 

C Ori A 

s 

XV 

0.77 

± 

1.21 

< 

1.66 

17 


— 



— 

7 

a Ori 

Si 

XIII 

0.21 

± 

0.37 

< 

1.02 

25 


— 



— 

17 

HD 206267 

Si 

XIII 

0.34 

± 

0.32 

< 

1.52 

59 

1.97 

± 

1.71 

> 

1.02 

25 

/ Per 

Si 

XIII 

2.35 

± 

1.06 

> 

1.90 

209 

0.30 

± 

0.26 

< 

1.12 

89 


“Measured 7 Z = f/i ratio with 68% confidence limits. 

^Inferred radius of formation (-/?*). 

“Total number of counts in the MEG and HEG spectra in the fit range considered 
(without rebinning or subtracting background). Note that this includes all three 
lines in the complex (r, i, and /), but that only about half of the counts are in the 
/ and i lines, which give the 7 Z ratio. 
d From Waldron et al. (2004). 


This table reports measured Tv ratios and inferred radii of formation from the erratum 
of WC07. These values are significantly different than those in the original article 


could vary 7Z about the best-fit value while allowing the 
other parameters to also vary to minimize the fit statistic. 
To obtain confidence intervals of 68.3%, 90%, and 95.4% one 
would find the interval on which the fit statistic increased 
by 1.0, 2.706, and 4.0, respectively. 

If one were using the normalizations of / and i as 
free parameters rather than 7Z, one could construct a two- 
parameter confidence region. Here the contours of 68.3%, 
90%, and 95.4% confidence would be defined by delta-fit- 
statistics of 2.30, 4.61, 6.17, respectively. Then one would 
find the values of the f/i ratio allowed from this plot. 

Because the confidence intervals reported by WC07 for 
7Z are symmetric about the best-fit value, it seems likely that 
they were not obtained using a delta-fit-statistic criterion, 
which will produce quite asymmetric confidence intervals 
(see § 4.5 for further discussion). Rather, they may have been 
computed using some other algorithm which finds a single 
relative error for each model parameter, perhaps based on 
the signal-to-noise ratio of the data. 

If such an algorithm was used, it is difficult to see how 
it could properly account for the covariance of the line nor- 
malizations with other parameters. Specifically, the derived 
normalization of the intercombination line is very sensitive 
to the strength and width assumed for the forbidden line. 
In order to properly account for this, it is necessary to make 
a two-parameter confidence region for the normalizations of 
the forbidden and intercombination lines while allowing all 
the other fit parameters to vary as well. 

A further potential problem with the methodology of 
WC07 is the use of the Doppler shift of the complex as 
a free parameter. Because the number of counts is small 
in many of the complexes considered, and the complexes 
have significant blending, a significant Doppler shift in the 
centroid can be degenerate with other model parameters. 
Furthermore, the centroid Doppler shift is expected to be 
close to zero, based on the profiles of stronger lines in the 
same part of the spectrum. I discuss this further in § 4.4. 


4 DATA REDUCTION AND ANALYSIS 

In § 4.1 I describe the reduction of the Chandra data used in 
this paper. In § 4.2 I describe the Gaussian fitting procedure 
used throughout this section. In § 4.3 I use this procedure 
to fit the data using a standard, robust procedure and find 
the 90% confidence intervals for all interesting parameters, 
as well as 68%, 90%, and 95% confidence intervals for 7Z. 
In the remainder of this section I discuss possible fitting 
artifacts and systematic errors. 

4.1 Data reduction 

All of the data considered in this paper are from archival 
Chandra HETGS (Canizares et al. 2000) observations of O 
stars. I give a list of all observations used in Table 3. 

The level one event Hies were reprocessed with 
acis_process_events to add up-to-date CTI corrections, af- 
ter which the data were reprocessed with standard CIAO 
routines for HETG spectra, producing spectrum and back- 
ground histograms, as well as response matrices and ancil- 
lary response files. CIAO 3.4 and CalDB 3.3.0 were used in 
all cases. 

Fitting was performed simultaneously on both gratings 
(HEG and MEG) on both positive and negative first orders 
and using all available observations, with no coadding or 
rebinning; the one exception to this is 6 1 Ori C, where I fit 
the two observations separately, as the magnitude of its X- 
ray emission is known to be rotationally modulated (Gagne 
et al. 1997). The use of unrebinned data is appropriate, since 
I use the C statistic and not y 2 (Cash 1979). 

The procedures described in the preceding paragraphs 
are essentially identical to those described in the data anal- 
ysis threads on the CIAO webpage 1 . 

1 http:/ /cxc. harvard.edu/ciao/threads/ 
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Table 3. List of observations with net exposure times 


Star 

obsid“ 

t b 

Cyg OB2 8A 

2572 

65.1 

C Pup 

640 

67.6 

e 1 Ori C 

3, 4 

49.5, 30.9 

C Ori A 

610, 1524 

59.6, 13.8 

a Ori 

3738 

91.0 

HD 206267 

1888, 1889 

34.1, 39.5 

£ Per 

4512 

158.8 


° Chandra Observation ID. 

^Net exposure time in ks. 

4.2 Gaussian fitting procedure 

In all cases, I used a model with five free parameters. Three 
parameters give the line normalizations: 7 Z = f /i, G = (f + 
i)/r, and the overall normalization of the complex. The other 
two parameters are the Doppler widths and shifts of the 
Gaussians, which were constrained to be the same for all 
lines in a complex. In all cases except where explicitly stated, 
the Doppler shift was fixed to zero. This is done to avoid 
fitting degeneracies, and because the shift is expected to be 
very close to zero for the complexes considered. I discuss the 
assumption of zero Doppler shift in more detail in § 4.4. 

The fitting procedure I used is similar to the approach 
taken by WC07, but with the important difference that I 
used the line ratios and overall complex normalization as free 
parameters rather than the individual line normalizations 
of /, i, and r. This does not significantly affect the best-fit 
values, but it greatly simplifies the extraction of confidence 
regions for the parameter of interest, TZ, using the delta-fit- 
statistic criterion (Lampton et al. 1976; Press et al. 2007). I 
emphasize that this procedure is identical to that of WC07 
in the sense that I am fitting three Gaussian emission lines 
with the same velocity broadening and shift, and the only 
difference is in the set of parameters I use as independent 
variables. 

Because / and i are blended due to both the finite reso- 
lution of the spectrometer and Doppler broadening, the nor- 
malizations of / and i have significant covariance. One will 
underestimate the uncertainty on TZ if one simply adds errors 
on the individual line normalizations in quadrature while ig- 
noring the covariance. However, when using TZ explicitly as a 
model parameter, its uncertainty can be computed directly 
by calculating the change in fit statistic obtained by varying 
it, thus bypassing the issue of evaluating the covariance of 
the individual lines. 

Another possible approach (not taken in this paper) 
would be to use the individual line normalizations as free 
parameters rather than the TZ and G ratios, and then to 
construct two-parameter confidence regions for the normal- 
izations of the forbidden and intercombination lines. The 
line ratio allowed for a given confidence level could be found 
from the allowed region on this plot. 

I do not subtract background in any fits. The most rig- 
orous way to include the effects of background would be 
to model the observed background spectrum and then add 
the modeled background into the total model for the source 
region. However, the number of background counts in the 
wavelength regions containing the He-like triplets is very 


small and can be adequately modeled as an addition to the 
weak bremsstrahlung continuum. As a representative exam- 
ple, based on the number of counts observed in the back- 
ground spatial region from 4-6 A, I estimate the number of 
background counts in the source extraction region of C, Ori 
near the S XV triplet to be 0.48 ±0.07 for the coadded first- 
order MEG data and 0.68 ± 0.08 for the coadded first-order 
HEG data, which is small compared with the observed num- 
ber of source counts in S XV of 17 in the MEG and 7 in the 
HEG. The other data sets have comparable or greater source 
count rates. Thus, I take the approach of treating the back- 
ground as a weak contaminant to the continuum, which is 
dominated by bremsstrahlung. 

I approximate the sum of the continuum and back- 
ground as a power-law of index two (which is flat when plot- 
ted in flux per unit wavelength vs. wavelength), with the nor- 
malization determined from a separate fit to the spectrum 
near the complex of interest; moderately strong lines are 
excluded from the continuum fit. This procedure is slightly 
different from that of WC07, who subtract the background 
and model the continuum as bremsstrahlung from a 10 MK 
plasma. However, given the small number of counts observed 
in each complex and the small spectral regions used in the 
fits, a power law of index two will characterize the contin- 
uum and background accurately enough to give equivalent 
results. 

In all cases, the He-like triplets were fit over a fixed 
range: 3.9-4.05 A for Ar XVII, 4.99-5.15 A for S XV, and 
6.59-6.79 A for Si XIII. 

4.3 Fit results 

I have fit all of the complexes given in Table 2 according to 
the procedure outlined in section 4.2. In Tables 4 and 5 I 
present the results of the fits. Table 4 gives best-fit models 
with 90% confidence limits for individual parameters. In Ta- 
ble 5 I present 68%, 90%, and 95% confidence intervals for 
TZ. 

The results presented in Tables 4 and 5 are in some sense 
a different measurement than that reported in WC07; I have 
performed a joint fit to the combined MEG and HEG data 
(but without coadding), while WC07 performed separate fits 
to the MEG and HEG data. In order to facilitate a direct 
comparison of results, I have also fit the MEG and HEG 
data separately. I report the best-fit values and the 68%, 
90%, and 95% confidence intervals on TZ in Tables 6 and 7. 
The results presented in Tables 4, 5, 6, and 7 are compared 
to the results of WC07 in Figures 1, 2, and 3. 

It should be noted that there is in general no astrophys- 
ical reason or instrumental to expect the MEG and HEG fits 
to give different results from each other or from a joint fit, 
and that I perform fits to the individual MEG and HEG 
data only in order to ensure a fair comparison of fit results. 
To the extent that one might expect instrumental differences 
between results from the different gratings, the best prac- 
tice is to identify specific cases where discrepancies exist, not 
just between the MEG and HEG, but between positive and 
negative dispersion orders. In cases where there is a discrep- 
ancy, one can attempt to resolve it or to discard data that 
may be compromised. In cases where there is no discrepancy, 
one should perform a joint fit to maximize the statistical ro- 
bustness of the results. In § 4.7 I examine individual grating 
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orders for the few cases that do show a discrepancy between 
HEG and MEG results. 

Although TZ is constrained to be less than TZo at 68% 
confidence from the joint fits in a few cases, I find that this is 
true for only one case at 90% confidence: the Si XIII complex 
of £ Per. I discuss this in more detail in § 4.8. Because all 
of the data considered (excepting £ Per) are consistent with 
1Z = TZo at 90% confidence, this implies that these data 
cannot provide an upper limit to the plasma radial location 
at 90% confidence. 

Comparing the measurements of TZ ratios made using 
only the MEG or only the HEG, I find significant, system- 
atic disagreement with the best-fit values and confidence 
intervals reported by WC07 for several cases. However, in 
two cases (the HEG measurement of Si XIII for £ Per, and 
the MEG measurement of Si XIII for a Ori) I find nominal 
agreement with the measurement by WC07 of a small value 
of TZ. These cases are considered in more detail in § 4.7. 

As a visual check on the validity of my models, I present 
the data for two He-like triplets as representative examples 
in Figures 4, 5, 6, and 7, as well as the best-6t models and 
models with TZ taking the extreme values allowed at 90% 
confidence, or AC = 2.706. It is difficult to believe that any 
of the models presented are excluded at 90% confidence or 
greater, as would be predicted by the reported confidence in- 
tervals of WC07. In particular, the model for the S XV triplet 
of Cyg OB2 8A with TZ = oo shows that it is possible to fit 
a model to this data set which has no intercombination line, 
with the broadening of the forbidden and resonance lines ac- 
counting for the observed counts near the wavelength of the 
intercombination line. According to the fit results of WC07, 
this model should be strongly excluded. (This should not 
be taken to imply that I believe that there may be no in- 
tercombination line flux, which would be a surprising result 
not predicted by atomic theory, but rather that the data are 
of poor statistical quality and do not even exclude a model 
with no intercombination line flux.) 


MEG/HEG 68% 

— MEG/HEG 90% 

MEG 68% 

HEG 68% 

WC07 MEG 68% 

WC07 HEG 68% 


Cyg OB2 8A 


<Pup 


COriA 


Figure 1. Comparison of several measurements of TZ for S XV 
observed in several stars. Measurements are at 68% confidence, 
with 90% confidence also given in gray for my joint fit to the HEG 
and MEG data. My measurements were made using the maximum 
likelihood method (C statistic). The dashed line indicates the 
value of TZo f° r S XV. The black points give TZ from my joint fits 
to MEG and HEG data, the red and magenta points give TZ from 
individual fits to MEG and HEG data, respectively, and the blue 
and cyan points give TZ from the individual fits of WC07 to MEG 
and HEG, respectively. 


5 


2 


MEG/HEG 68% 

MEG/HEG 90% 
MEG 68% 

— HEG 68% 

WC07 MEG 68% 

WC07 HEG 68% 


4.4 The effects of assumptions on velocity 
broadening and shifts on fit results 

The fact that in some cases I have measured significantly 
different best-fit values and confidence intervals for TZ than 
WC07 demands explanation. Unfortunately, WC07 have not 
reported enough information to quantitatively reproduce 
their results. Although they report best-fit values with 68% 
confidence limits for TZ = f/i and G = (/ + i)/r, they do 
not report velocity widths or shifts. However, it is possible 
to test hypotheses regarding the origin of some of the dis- 
crepancies. 

The greatest discrepancy between my measurements 
and those of WC07 is in the S XV triplet of Cyg OB2 8A. A 
clue to a possible source for this disagreement can be found 
in their reported G = (/ + i)/r ratios. Using the MEG, 
WC07 measured G = 3.16 ± 1.70, and using the HEG, they 
measured G = 1.98± 1.57 (both at 68% confidence). In com- 
parison, I measured G = 0.81q3 at 90% confidence. The 
best-fit value I found is much closer to the range expected 
from atomic physics ( G ~ 1.0, e.g. Pradhan 1982; Porquet 
et al. 2001). Furthermore, the visual appearance of the data 
(shown in Figs. 4 and 5) suggests that G ~ 1 and not that 
G ~ 3. 


u Ori HD 206267 ( Per 

Figure 2. Same as Figure 1, but for Si XIII. 

One possible way to get a larger value for G would be 
to assume a significant blueshift in the line wavelengths, in 
which case one would be effectively counting photons near 
the rest wavelength of the resonance line as intercombination 
line photons. This would then increase the model value of i 
and decrease the model value of r. As discussed in § 4.2, I 
have assumed zero velocity shift in my fits. I have done this 
because one expects velocity shifts of zero for lines where 
the wind optical depth is small, and because I wanted to 
minimize the number of free parameters to avoid degeneracy. 

To test this idea, I have fit models to the MEG and 
HEG data where TZ and G were forced to take on the values 
reported by WC07. I allowed both the velocity width and the 
velocity shift to be free. The results are reported in Table 8. 
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Table 4. Best fit models for He-like triplets using maximum likelihood method 


Star 

Ion 

77 

G 

o a 

n h 

c 

bins 

Cyg OB2 8 A 

Ar XVII 

( 0 . 8 ) 

0 +ae 

950+ 850 

you _400 

8 + 2.8 

159.5 

176 

Cyg OB2 8 A 

S XV 

> 0.8 (2.4) 

o-stgj 

950+ 250 

you -200 

41±« 

196.5 

188 

C Pup 

S XV 

0 9 +1 ' 5 

U ‘ J -0.5 

0-8to.3 

700 +20 ° 
' uu — 200 

26±| 

174.8 

188 

9 1 Ori C (obsid 3) 

S XV 

i O+0.8 

— 0.4 

0 - 8 to .2 

350~*~ 100 

oou -100 

99±i° 

265.4 

188 

9 1 Ori C (obsid 4) 

S XV 

1 4 +1 ' 4 

0 - 8 to .3 

299tjjg 

76tie 

205.2 

188 

C Ori A 

S XV 

< 4.5 (0.2) 

0 4 + 0 - 8 
u -^-0.3 

Q50+ 350 

OOU -250 

2-5i{1 

111.7 

376 

a Ori 

Si XIII 

1 1 +1 - 7 
± - 1 - 0.7 

1 O " 1 ” 0 ' 9 
i - U - 0.5 

< 350 (200) 

9 q + 0.9 
^'- 0.6 

119.7 

236 

HD 206267 

Si XIII 

> 0.5 (1.4) 

Q ' 7 + 0.5 

u ‘ ' — 0.4 

700~*~ 300 

' uu — 250 

5-5±i;i 

256.0 

472 

§ Per 

Si XIII 

1 l +0 - 7 
± - i - 0.4 

i O+0.5 
0.3 

700" 1 ” 200 

' uu -100 

9.5±i;i 

263.6 

236 


a Best-fit Gaussian velocity width in units of kms . 
b Best-fit normalization in units of 10 -6 photons cm -2 s — 1 . 


All confidence intervals reported in this table are 90%, or AC = 2.706. Parenthe- 
ses denote a value that minimizes the fit statistic C where the parameter is either 
unconstrained at 90% confidence or gives only an upper or lower limit. 


Table 5. Allowed ranges for 7 Z at different confidence levels 


Star 

Ion 

best fit 

68.3% 
min max 

90% 

min max 

95.4% 
min max 

Cyg OB2 8 A 

Ar 

XVII 

0.8 

0.1 

3.3 

0 

OO 

0 

OO 

Cyg OB2 8 A 

S 

XV 

2.4 

1.2 

10.7 

0.8 

oo 

0.7 

oo 

C Pup 

S 

XV 

0.9 

0.5 

1.5 

0.4 

2.4 

0.3 

3.5 

9 1 Ori C (obsid 3) 

s 

XV 

1.2 

0.9 

1.6 

0.8 

2.0 

0.7 

2.2 

9 1 Ori C (obsid 4) 

s 

XV 

1.4 

1.0 

2.2 

0.8 

2.9 

0.7 

3.4 

C Ori A 

s 

XV 

0.2 

0 

1.0 

0 

4.5 

0 

OO 

cr Ori 

Si 

XIII 

1.1 

0.6 

1.9 

0.4 

2.8 

0.3 

3.7 

HD 206267 

Si 

XIII 

1.4 

0.7 

3.2 

0.5 

OO 

0.4 

OO 

£ Per 

Si 

XIII 

1.1 

0.9 

1.5 

0.7 

1.8 

0.7 

2.0 


5 


3 


Pi 

2 


0 


MEG/HEG 68% 

— MEG/HEG 90% 

MEG 68% 

— HEG 68% 

WC07 MEG 68% 

WC07 HEG 68% 


i 


obsid 3 


obsid 4 


both 


Figure 3. Same as Figure 1, but showing results for S XV in 
the spectrum of 9 1 Ori C. My fits were done for obsid 3 and 4 
separately, while those of WC07 were done for coadded data. 


I find that models requiring their reported best-fit values of 
77 and G indeed prefer large blueshifts. 

Since many emission line profiles in the X-ray spectra 
of O-stars are indeed quite blueshifted, especially for stars 


with high mass-loss rates, it is reasonable to ask whether the 
S XV triplet of Cyg OB2 8A might in fact be significantly 
blueshifted. One would expect an insignificant blueshift be- 
cause of the relatively small atomic opacity of the wind at 
5 A. However, the most reliable way to test the expected 
blueshift is to fit the profile of a strong line nearby in the 
spectrum. 

In Table 9, I report the results of a Gaussian fit to 
the Si XIV Ly a line of Cyg OB2 8A. I indeed find that 
the blueshift is consistent with zero, and that the velocity 
broadening is consistent with the measured value from my 
fit to S XV. Therefore, at least in this case, my assumptions 
about the velocity structure of S XV appear to be justified, 
and any fit to S XV which prefers a large blueshift must be 
regarded as an artifact. 

Although it is not possible to know what fitting pro- 
cedure led WC07 to obtain the results they report for Cyg 
OB2 8A, it seems reasonable to assume that allowing the ve- 
locity shift to be free may have led to degeneracy in model 
parameters. However, I have shown that fixing the velocity 
shift to zero is a good assumption, and that under this as- 
sumption, I find much more reasonable best-fit values for G. 
If this is indeed the case, then the low values of 77 reported 
by WC07 for S XV in the spectrum of Cyg OB2 8A are likely 
an artifact of a fitting degeneracy. 
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Table 6. Allowed ranges for 1Z from fitting only the MEG data 


Star 

Ion 

best fit 

68.3% 
min max 

90% 

min max 

95.4% 
min max 

Cyg OB2 8A 

Ar 

XVII 

0.2 

0 

1.0 

0 

OO 

0 

OO 

Cyg OB2 8A 

S 

XV 

4.5 

1.6 

OO 

1.0 

OO 

0.8 

OO 

C Pup 

S 

XV 

0.9 

0.5 

2.2 

0.3 

8 

0.2 

OO 

e 1 Ori C (obsid 3) 

s 

XV 

1.3 

0.9 

1.9 

0.7 

2.6 

0.7 

3.2 

d 1 Ori C (obsid 4) 

s 

XV 

3.5 

2.1 

6.6 

1.6 

12 

1.3 

23 

C Ori A 

s 

XV 

0.4 

0.1 

1.4 

0 

3.5 

0 

6.2 

<j Ori 

Si 

XIII 

0.4 

0.2 

0.8 

0.1 

1.3 

0.1 

1.8 

HD 206267 

Si 

XIII 

1.0 

0.5 

3.4 

0.3 

OO 

0.2 

OO 

§ Per 

Si 

XIII 

1.9 

1.4 

2.9 

1.1 

4.2 

1.0 

5.6 


Table 7. Allowed ranges for 7 Z from fitting only the HEG data 


Star 

Ion 

best fit 

68.3% 
min max 

90% 

min max 

95.4% 
min max 

Cyg OB2 8A 

Ar XVII 

1.9 

0.2 

OO 

0 

OO 

0 

OO 

Cyg OB2 8A 

S XV 

1.3 

0.5 

4.8 

0.2 

OO 

0 

OO 

C Pup 

S XV 

1.1 

0.6 

2.2 

0.4 

6.0 

0.2 

OO 

d 1 Ori C (obsid 3) 

S XV 

3.0 

1.3 

OO 

0.9 

OO 

0.8 

OO 

6 1 Ori C (obsid 4) 

S XV 

0.6 

0.2 

1.1 

0.1 

1.7 

0.1 

2.0 

C Ori A 

S XV 

0 

0 

OO 

0 

OO 

0 

OO 

<j Ori 

Si XIII 

6.7 

2.1 

OO 

1.2 

OO 

0.9 

OO 

HD 206267 

Si XIII 

1.5 

0.8 

3.5 

0.5 

7.0 

0.4 

12 

£ Per 

Si XIII 

0.3 

0.1 

0.7 

0.0 

1.0 

0.0 

1.3 



Wavelength (A) 



Wavelength (A) 


Figure 4. First order MEG spectrum of Cyg OB2 8A (black hor- 
izontal lines with error bars) shown near the wavelengths of the 
S XV triplet together with three models. The positive and nega- 
tive first orders have been coadded and the spectrum rebinned by 
a factor of two for presentation purposes. The rest wavelengths 
of the lines are indicated with dashed vertical lines. The error 
bars on the individual bins are for presentation purposes and are 
calculated using the prescription of Churazov et al. (1996). Note 
that the C statistic (used in this paper) is not computed using a 
variance function, unlike x 2 (Cash 1979). The models are from a 
joint fit with HEG data. The best-fit model is in red, and models 
with lower and upper limits to the TZ parameter at 68% confidence 
are in blue and green, respectively. 


Figure 5. As figure 4, but showing HEG data for S XV in the 
spectrum of Cyg OB2 8A with models. The positive and negative 
first orders have been coadded, and the spectrum rebinned by a 
factor of four for presentation purposes. 


4.5 Asymmetry of confidence limits on 7 Z for 
S XV in the spectrum of £ Pup 

I find moderate disagreement with WC07 in the size of the 
error bars and the best-fit values of 7Z obtained using fits to 
only the MEG or HEG data. I also find that the confidence 
regions on 7 Z are quite asymmetric, in contrast with the 
symmetric confidence regions reported by WC07. This is 
true for all the data sets, but it is especially important for 
the S XV triplet in the spectrum of £ Pup. 
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Wavelength (A) 


Figure 6. As figure 4, but showing MEG data for S XV in the 
spectrum of £ Pup with models. The positive and negative first 
orders have been coadded, and the spectrum rebinned by a factor 
of two for presentation purposes. 



Wavelength (A) 


Figure 7. As figure 4, but showing HEG data for S XV in the 
spectrum of £ Pup with models. The positive and negative first 
orders have been coadded, and the spectrum rebinned by a factor 
of four for presentation purposes. 


Table 8. Best-fit velocity parameters assuming line ratios re- 
ported by WC07 for S XV in the spectrum of Cyg OB2 8A 


Grating 

TZ 

G 

cr a 

Av b 

n“ 

MEG 

(1.07) 

(3.16) 

1450 

-725 

39 

HEG 

(1.13) 

(1.98) 

325 

-700 

41 


“Best- fit Gaussian velocity width in units of kms -1 . 
^Best-fit velocity shift in units of kms -1 . 

“Best- fit normalization in units of 10 — 6 photons cnr 2 s -1 . 

Parentheses denote fixed parameters taken from WC07. 


Table 9. Best-fit Gaussian model for Si XIV in the spectrum of 
Cyg OB2 8A 


cr a 

•O 

<1 

n“ 

850ll® 0 

-50 ±100 

34 ±4 


“Best-fit Gaussian velocity width in units of kms -1 . 

Best.- fit velocity shift in units of kms -1 . 

“Best-fit normalization in units of 10 -6 photons cm -2 s -1 . 

Fit results are for a joint fit to MEG and HEG data. Confidence 
intervals are reported at 90%. 

The asymmetry is not surprising, as the dependence of 
the fit statistic on the line ratio is expected to be nonlinear. 
As the line ratio becomes very large or very small, the data 
are no longer very sensitive to changes in the model ratio. 

As discussed in § 3.2, it is not clear from the text of 
the erratum of WC07 how the error bars they reported were 
derived. However, the symmetry of the error bars reported 
by WC07 suggests that they may have used a simplified or 
erroneous procedure for estimation of parameter confidence 
intervals. It also implies that their procedure did not use a 
delta-fit-statistic criterion to evaluate parameter confidence 
intervals. 

4.6 The Si K edge in the ACIS detector 

The effective area of the front-illuminated chips of the ACIS- 
S detector used in reading out the Chandra HETGS shows 
a strong Si K edge which has significant calibration uncer- 
tainties. There is also a Si K edge in the effective area of the 
back-illuminated chips, but it is much less pronounced. 

When using the standard aimpoint for the location of 
a bright source in the focal plane, the wavelength of the Si 
K edge falls on a back-illuminated chip for the positive first 
order MEG data, while it falls on a front-illuminated chip for 
the negative first order MEG data, as well as both positive 
and negative first order HEG data. This results in a strong Si 
K edge in the effective area of all first order HETGS spectra 
except the MEG positive first order. This is illustrated in 
Fig 8, where I plot the effective area of all four first order 
grating arms in the HETGS spectrum of £ Per in the vicinity 
of the Si K edge. 

The fact that the positive first-order MEG data do not 
have a strong Si K edge implies that they are the best data 
for measuring 1Z for Si XIII, and that if there is any dis- 
crepancy between any of the grating orders, the positive 
first-order MEG data should be considered the most reli- 
able. This is in contradiction to the claim of WC07 that the 
HEG data have a weaker Si K edge than the MEG data. 

4.7 Fits to individual orders: Si XIII in a Ori and 
£ Per 

I have measured quite low values for 1Z in two cases: using 
only the MEG for Si XIII in the spectrum of a Ori and using 
only the HEG for Si XIII in the spectrum of £ Per. In each 
of these cases, the joint fit to both the HEG and MEG data 
gives a much higher value for 7 Z. Both of these measurements 
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— ±1 



+1 68% 



-1 68% 



WC07 ±1 68% 
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1 


a Ori MEG ( Per HEG 


Figure 8. Effective area of individual gratings and orders of 
HETGS near the wavelengths of the Si XIII triplet for the ob- 
servation of £ Per. The rest wavelengths of the triplet are shown 
with vertical dashed lines. There is a prominent Si K edge near 
the wavelength of the forbidden line of Si XIII for all the orders 
which are read out on front-illuminated chips (negative first order 
MEG, and both first order HEG), but only a small edge for the 
positive first order MEG, which is read out on a back-illuminated 
chip. 


to individual gratings agree with the results reported by 
WC07. Both also disagree at 68% confidence with the results 
obtained from the same observation, but using the other 
grating (HEG instead of MEG, or vice versa), although they 
are formally consistent at 95% confidence. The disagreement 
between different gratings calls for closer examination. 

In order to assess the validity of these measurements I 
have fit individual grating arms (positive and negative first 
orders). Although there is in general no reason to expect 
different results from the individual orders, there are two 
reasons why this might be useful: first, when looking at pos- 
itive and negative first order MEG data, uncertainties in the 
calibration of the strong Si K edge feature in the negative 
first order MEG data (described in § 4.6) may lead to dis- 
crepancies in the two orders; and second, fits to individual 
orders serve as a consistency check on the results from the 
joint fit, helping to find artifacts or statistical flukes. 

In order to avoid scrutinizing only data for which I have 
perceived an unexpected outcome, I have performed these 
fits for all data sets considered in this paper. However, I only 
found significant disagreement between positive and nega- 
tive first orders for these two cases. I report these fit results 
in Table 10. The results are also depicted in comparison with 
those of WC07 in Figure 9. 

For the Si XIII data of a Ori, I find that although the 
joint fit to the MEG data indicates a quite low value of 7 7, 
the two individual grating orders show very different values, 
although they are formally consistent at the 90% confidence 
level. The negative first order indicates that 77 is relatively 
small, while the positive first order indicates a much larger 
value. Since the negative first order may be affected by the 
strong Si K edge of the front-illuminated CCD chip it falls 
on, this marginal discrepancy may be resolved by consid- 
ering the positive first order data as the most reliable, as 


Figure 9. Same as Figure 2, but showing results for individual 
grating orders for the Si XIII complexes of two specific stars which 
show discrepancies between orders. 

discussed in § 4.6. Furthermore, the joint fit to the HEG 
data also indicates a higher value for 77. 

For the Si XIII data of £ Per, I found that the joint best 
fit to the HEG data required an unrealistically high veloc- 
ity broadening, such that the lines of the triplet were only 
marginally resolved from one another. Again, I find that 
the individual grating orders show very different values of 
77, but are formally consistent at the 90% confidence level. 
However, I also find that the fit to the negative first order 
HEG data, which prefers a higher value for 77, has a reason- 
able value for the velocity broadening, which is comparable 
to the broadening obtained by fitting the MEG data, or to 
the broadening obtained by fitting lower Z lines (with more 
counts) in the same spectrum . On the other hand, the fit 
to the positive first order HEG data, which prefers a low 
value for 77, requires an unrealistic value for the velocity 
broadening, as I found for the joint fit. Thus, while there is 
no reason to prefer one order over the other, it seems rea- 
sonable to suppose that the discrepancy between the two 
orders, combined with the strange velocity distribution of 
counts in the positive first order HEG data, indicates that 
these data are either a statistical fluke, or in some other 
way anomalous. Furthermore, the MEG data are consistent 
with the hypothesis that 77 = 77o- Since the MEG data are 
unambiguous, and the HEG data show a discrepancy, it is 
reasonable to conclude that the data are at least as consis- 
tent with the hypothesis that 77 = 77o as with any other 
hypothesis. 

4.8 Constraints on plasma formation radius from 
the joint fit to Si XIII in the spectrum of £ Per 

Although the constraints on the 77 ratio of Si XIII in the 
spectrum of £ Per may be skewed by the possibly anoma- 
lous positive first order HEG data, as discussed in § 4.7, 
I will use the formal measurement of 77 from the joint fit 
to MEG and HEG data to derive constraints on the plasma 
formation radius under the assumption of a single formation 
radius, in combination with the calculations of Leutenegger 
et al. (2006). I do this by comparing the measured 77 ratios 
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Table 10. Comparison of allowed ranges for 77 from fitting individual grating orders 
in two selected cases 


Star 

Ion 

Grating 

best fit 

68.3% 

min max 

90% 

min max 

95.4% 
min max 

o 

b 

Si XIII 

MEG ±1 

0.4 

0.2 

0.8 

0.1 

1.3 

0.1 

1.8 



MEG +1 

1.3 

0.6 

3.8 

0.3 

12 

0.2 

OO 



MEG -1 

0 

0 

0.1 

0 

0.8 

0 

1.1 

£ Per 

Si XIII 

HEG ±1 

0.3 

0.1 

0.7 

0.0 

1.0 

0.0 

1.3 



HEG +1 

0.2 

0.0 

0.3 

0.0 

0.4 

0.0 

0.6 



HEG -1 

2.3 

1.1 

OO 

0.1 

OO 

0 

OO 


from this work to the plots of 1Z(R) for £ Pup and l Ori in 
Figures 3 and 5, respectively, of Leutenegger et al. (2006). 
The photospheric temperature and spectral type of £ Per is 
intermediate between those of £ Pup and l Ori, so the be- 
havior of 1Z for £ Per should be bracketed by that of £ Pup 
and i Ori. 

Using the measured 90% confidence limits on the 1Z 
ratio of Si XIII in the spectrum of £ Per (0.7 < 1Z < 1.8), 
I find an implied plasma location of 2 < R < 6 using the 
£ Pup model UV flux, and 1.4 < R < 4 using the l Ori model 
UV flux. The ranges given here are statistical uncertainties, 
not plasma distribution ranges. In either case, the data do 
not require an X-ray emitting plasma located very close to 
the photosphere. The location inferred is comparable to that 
inferred for other O stars under the assumption of a single 
plasma formation radius. 

4.9 Comparison with previous results for £ Pup 

The S XV triplet of £ Pup has been studied in several pre- 
vious works. I compare the previous and current results in 
Figure 10. 

My current results are consistent with the measure- 
ments of Leutenegger et al. (2006). It is not possible to 
comment on the consistency of the results of WC07 and 
their earlier work, since Cassinelli ct al. (2001) did not re- 
port which gratings or what confidence intervals were used. 
Again, while 7Z = 77o is ruled out by my measurements at 
the 68% confidence level, it is allowed at 90%. Thus, these 
data provide no upper limit on the radial location of the 
X-ray emitting plasma. 

Even using the 68% confidence upper limit of 77 < 1.5, I 
only infer an upper limit to the radius of S XV of 2.2R» using 
Figure 3 of Leutenegger et al. (2006) under the assumption 
of a single plasma formation radius. This limit is completely 
consistent with the numerical modeling results dicussed in 
the introduction to this paper. 


5 DISCUSSION 

In § 4.3 I have measured the 7 Z = f /i ratios of high-Z He- like 
triplets of several O stars with a simple, robust procedure. 
I find that in all but one of the cases considered, 77 is con- 
sistent with 77o at 90% confidence. The one case where 77 is 
inconsistent with 77o (Si XIII on £ Per) is expected, based 
on the strong photospheric UV flux of that star, and does 
not imply that the formation radius must be especially close 
to the photosphere. 



68% confidence 

90% confidence 

LPKC06 68% 

WC07 MEG 68% 

— WC07 HEG 68% 
C+01 (68% ?) 




« 


0.5 ■ 

0.0 

C Pup S XV 

Figure 10. Comparison of results from this work with all previ- 
ous results on the same data for S XV in the spectrum of £ Pup. 
LPKC06 refers to Leutenegger et al. (2006), and C+01 refers to 
Cassinelli et al. (2001). 


The fact that I find that 77 is consistent with 77o for 
most of the cases considered leads directly to the conclu- 
sion that none of those measurements can provide an upper 
limit to the radius of formation, and cannot be used to in- 
fer the necessary presence of X-ray emitting plasma close to 
the stellar photosphere. Because the data considered in this 
paper are of poor statistical quality, the location of the X- 
ray emitting plasma is unconstrained or poorly constrained, 
and the possibility of X-ray emitting plasma near the photo- 
sphere is not excluded; however, the inference of X-ray emis- 
sion close to the photosphere would be a surprising result 
requiring strong positive evidence to establish credibility. 

Since there is a large body of evidence indicating that 
all of the X-ray emitting plasma other than the high-Z ions 
is distributed throughout the winds of O stars, including the 
velocity broadening of the emission lines and the f /i ratios 
of low-Z ions (e.g. Waldron & Cassinelli 2001; Kahn et al. 
2001; Cassinelli et al. 2001), it would be natural to infer that 
the high-Z ions are formed in shocks distributed throughout 
the winds as well, unless there is significant evidence to the 
contrary. In this work I have shown that the f /i ratios do not 
constrain these high-Z ions’ locations. However, WC07 have 
shown that in a number of cases high-Z ions inferred by them 
to be close to the photosphere showed significant velocity 
broadening, which would not be expected for lines forming 
very close to the photosphere, where the mean flow velocity 
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is small. I have confirmed their result: many of these lines 
do indeed show significant velocity broadening, as shown in 
Table 4. Although WC07 framed the discrepancy between 
their inferred plasma locations and the large velocity broad- 
enings of the lines as a new, interesting problem, the data 
are in fact completely consistent with the hypothesis that 
the higli-Z ions are not forming near the photosphere. 

The fact that the data considered in this paper do not 
significantly constrain the location of the high-temperature 
X-ray emitting plasma does not invalidate the conclusion 
that other He-like triplets do provide meaningful constraints 
(Kahn et al. 2001; Waldron & Cassinelli 2001; Cassinelli 
et al. 2001; Miller et al. 2002; Waldron et al. 2004; Leuteneg- 
ger et al. 2006, 2007). However, in all of those cases, lower-Z 
ions were found to be at locations in the wind that were 
consistent with a wind-shock origin for the X-ray emitting 
plasma. Even in the case of 9 1 Ori C, which has a mag- 
netically channeled wind, the X-ray emitting regions may 
be significantly removed from the photosphere (Gagne et al. 
2005). 


6 SUMMARY 

I have attempted to reproduce measurements of 1Z = f/i by 
WC07 for a number of stars which led them to infer plasma 
formation radii close to the stellar photospheres. I find that 
all but one of the complexes considered is consistent with 
1Z = IZo at 90% confidence, and thus with no photoexcita- 
tion effect; the one complex that has 1Z < IZo is expected 
based on the strong UV flux of that star, and does not re- 
quire plasma formation near the photosphere. Because the 
f/i ratios examined in this paper are all poorly constrained 
or unconstrained by the data, in most cases the data provide 
no interesting constraint on the location of the X-ray emit- 
ting plasma. There is no evidence requiring the existence 
of X-ray emitting plasma very close to the photospheres of 
O stars, and there is no evidence for a near-star high-ion 
problem, as reported by WC07. 

I have made these measurements using a robust proce- 
dure with an orthodox determination of the 68%, 90%, and 
95% confidence intervals. I used a delta-fit-statistic criterion 
to compute these confidence intervals, so that model param- 
eter covariance is properly taken into account. In particular, 
the high-Z He-like triplets are significantly blended, leading 
to non-negligible covariance in the line normalizations. 

I have also shown that, for data of marginal statistical 
quality, if the centroid shift of the complex is taken as a free 
parameter, it may be degenerate with the line normaliza- 
tions. By measuring the centroid shifts of nearby emission 
lines of good statistical quality, I have shown that it is jus- 
tified to assume zero centroid shift for the high-Z he-like 
complexes and remove this degeneracy. 

These results on data of low statistical quality from 
high-Z He-like ions do not contradict or invalidate the many 
earlier results finding meaningful constraints on the loca- 
tion of lower-Z ions using data of higher statistical quality. 
The low-Z ions are located in the winds of O stars, and the 
high-Z ions’ locations are not constrained by the data. 

WC07 have measured significant velocity broadening 
for a number of high-Z ions, and I confirm this result. In 
the light of the lack of constraint of the radial location of 


these ions by their f/i ratios, this implies that these ions 
are formed in the wind and not near the photosphere, as is 
the case for the rest of the X-ray emitting plasma. 
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